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We comment on the results of a complete leading-twist next-to-leading order QCD 
analysis of the spacelike pion electromagnetic form factor at large-momentum 
transfer Q. For the asymptotic distribution amplitude, we have examined the 
sensitivity of the predictions to the choice of the renormalization scale. The re- 
sults show that reliable perturbative predictions for the pion electromagnetic form 
factor can already be made at a momentum transfer Q below 10 GeV. 
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1. Introduction 

The framework for analyzing exclusive processes at large-momentum transfer within 
the context of perturbative QCD (pQCD) has been developed in the late seven- 
ties Q. It was demonstrated, to all orders in perturbation theory, that exclusive 
amplitudes involving large-momentum transfer factorize into a convolution of a 
process- independent and perturbatively incalculable distribution amplitude (DA), 
one for each hadron involved in the amplitude, with a process-dependent and per- 
turbatively calculable hard-scattering amplitude. In the standard hard-scattering 
approach (sHSA), hadron is regarded as consisting only of valence Fock states, 
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transverse quark momenta are neglected (collinear approximation) as well as quark 
masses. 

Although the pQCD approach of Ref. undoubtedly represents an adequate 
and efficient tool for analyzing exclusive processes at very large momentum trans- 
fer, its applicability to these processes at experimentally accessible momentum 
transfer has been long debated and attracted much attention. There are several 
important issues regarding this subject. Let us mention first that in a moderate 
energy region (a few GeV) soft contributions (resulting from the competing, so- 
called, Feynman mechanism) can still be substantial J2j, but the estimation of their 
size is model dependent. Further on, the self-consistency of the pQCD approach 
was questioned regarding the nonfactorizing end-point contributions ||. It has 
been shown, that the incorporation of the Sudakov suppression in the so-called 
modified perturbative approach (mHSA) j|, ^ effectively eliminates soft contribu- 
tions from the end-point regions and that the pQCD approach to the pion form 
factor begins to be self-consistent for a momentum transfer of about Q 2 > 4 GeV 2 
[Q. However, in the pQCD approach to exclusive processes one still has to check 
its self-consistency by studying radiative corrections. 

It is well known that, unlike in QED, the leading-order (LO) predictions in 
pQCD do not have much predictive power, and that higher-order corrections are 
important. In general, they have a stabilizing effect reducing the dependence of 
the predictions on the schemes and scales. Although the LO predictions within the 
sHSA (as well as, the mHSA) have been obtained for many exclusive processes, 
only a few processes have been analyzed at the next-to-leading order (NLO): the 
pion electromagnetic form factor || ^, |), the pion transition form factor (7[ ||] 
(and Q in mHSA), and the process 77 -> MM (M = tt, A") (T^j. 

In our recent work Q we have clarified some discrepancies between previous 
results ||, and by including the complete closed form for the NLO evolution of 
the pion DA derived recently , we have obtained the complete NLO prediction 
for the pion electromagnetic form factor. 

In this work we would like to give a short summary of our calculation (for 
details and notation see ^|) with special emphasis on the proper choice of the 
renormalization scale. 

2. Pion electromagnetic form factor in the sHSA 

In leading twist, the pion electromagnetic form factor can be written as 

F 7I (Q 2 )= f dx [ dy $*(y,v 2 F )T H (x,y,Q 2 ,[i 2 R ,fi 2 F ) ${x,nl). (1) 
Jo Jo 

Here Q 2 = —q 2 is the momentum transfer in the process and is supposed to be 
large, fj,R is the renormalization scale, and /if is the factorization scale at which 
soft and hard physics factorize. 
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The hard-scattering amplitude Th{x, y, Q 2 , /ifj, fi F ) is the amplitude for a par- 
allel q{q 2 P a i r of the total momentum P, with the constituents carrying the mo- 
mentum fractions x and x = 1 — x, hit by a virtual photon 7* of momentum q 
to end up as a parallel qiq 2 pair of momentum P' = P + q, with the constituents 
sharing fractions y and y = 1 — y. It can be calculated in perturbation theory and 
represented as a series in the QCD running coupling constant a,s(/ifj): 



Tnix^.Q 2 , H R ,(J. F ) = a s{^ R )T^'{x,y,Q 



T^ > {x 1 y 1 ix R /Q z 1 ^ F /Q z ) + 



(2) 



There are, up to the order we are calculating, 4 LO and 62 NLO Feynman dia- 
grams that contribute to the (9192) + 7* ~~ * (Sifc) amplitude. We have used the 
dimensional regularization method and the MS renormalization scheme in our 
calculation (details and results are presented in B). 

The intuitive interpretation of the pion DA $(x,fj, R ) /ufj)) is that it rep- 

resents a probability amplitude for finding the valence qiq 2 Fock state in the initial 
(final) pion. The function <!> is intrinsically nonperturbative, but its evolution can 
be calculated perturbatively. It is advantageous to express the DA (further on, we 
use the function (f> normalized to unity) in terms of the Gegenbauer polynomials 



(x) = (p(x,nl) = (f> as {x) 



l + J2'B n C 3 J 2 (2x-l) 



n=2 



(3) 



Here cf) as (x) = 6x(l — x) is the asymptotic DA which represents the solution of 
the DA evolution equation for pr F — > 00, while the coefficients B n (n even) are 
obtained using some nonperturbative techniques at energy The DA is then 
evoluted to the energy \j? F and the DA evolution up to the NLO has the form 



(x,n%) 



(4) 



where we turn to @] for detailed expressions. 

Two most exploited choices for the pion DA (||) are <fi a s{x) and <ficz(x), for 
which E>2 = and B 2 = 2/3, respectively, while B n = for n > 2. Unlike tp as (x), 
the <pcz(x) is a strongly end-point concentrated distribution and its form has been 
obtained using the method of QCD sum rules . There is no LO evolution for 
4> as (x, (J, F ) and the NLO evolution is tiny. As we have shown in the inclusion 
of the LO evolution is crucial when one tends to obtain meaningful results with 
the 4>cz function, and even the NLO evolution is significant. 

Although the numerical results for the pion electromagnetic form factor ob- 
tained by using the 4>cz{x, fi F ) distribution are higher and closer to the existing 
experimental data, there are compelling theoretical results which disfavor the (j>cz 
distribution: theoretical predictions for the pion transition form factor are in very 
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good agreement with the data assuming that the pion distribution amplitude is 
close to the asymptotic one, i.e., 4> as {x, [M F ) |Q; the estimation of the size of the 
soft contributions to the pion electromagnetic form factor indicates the asymp- 
totic form of the pion DA, and even the self-consistency of the derivation of the 
4>cz(x, (ip) distribution from the QCD sum rules was criticized Taking this 
into account, one expects that the pion DA does not differ much from the <j> as and 
in this work we comment only on the results obtained with <p as (x, fip). 

By inserting (0) and (0) into ([!]) one obtains the complete NLO pQCD expres- 
sion for the pion electromagnetic form factor. Generally, one can express the NLO 
form factor as 

F„(Q 2 , v% /4) = FM (Q 2 , A) + (Q 2 , /£, Mf) • (5) 

The first term in (|J) is the LO contribution, while the second term is the NLO 
contribution coming from the NLO correction to the hard-scattering amplitude as 
well as arising from the inclusion of the NLO evolution of the DA. For the results 
obtained using (j> as (x, fi F ) distribution, the effect of the NLO evolution of the DA 
is negligible (w 1%). 



3. Choosing the factorization and renormalization scales 

The physical pion form factor Ft T (Q 2 ), represented at the sufficiently high Q 2 by 
the factorization formula (Q), is independent of the renormalization scheme and of 
the renormalization and factorization scales, hr and fj,p , respectively. Truncation 
of the perturbative series of F 7V (Q 2 ) at any finite order causes a residual dependence 
on the scheme as well as on the scales (which is already denoted in Eq. (^)). We 
approximate F^{Q 2 ) only by two terms of the perturbative series and hope that 
we can minimize higher-order corrections by a suitable choice of (1r and fj,p, so 
that the LO term F^{Q 2 , jjL 2 Rl ii 2 F ) gives a good approximation to the complete 
sum F„(Q 2 ). 

The simplest and widely used choice for the fin and \if scales is 

a4 = = Q 2 . ( 6 ) 

the justification for the use of which is mainly pragmatic. Physically, a more 
appropriate choice for \j? R would be that corresponding to the characteristic virtu- 
alities of the particles in the parton subprocess, which is considerably lower than 
the overall momentum transfer Q 2 (i.e., virtuality of the probing photon). The 
physically motivated choices we are using are 

n\ = xyQ 2 , (7) 

A = V(xyQ 2 )(VQ 2 ), (8) 

/4 = e^l^xyQ 2 . (9) 
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These correspond, respectively, to the (LO) gluon virtuality, geometrical mean of 
the gluon and quark virtualitics (an attempt to take into account that in QCD, 
unlike in QED, the coupling is renormalized not only by the vector particle prop- 
agator, but also by the quark-gluon vertex and the quark-propagator, Dittes and 
Radyushkin M ) , and to the choice of the renormalization scale according to the 
Brodsky-Lepage-Mackenzie (BLM) procedure ]l6| . The essence of the BLM pro- 
cedure is that all vacuum-polarization effects from the QCD (3 function are re- 
summed into the running coupling constant, and, as a result of the choice (^]), 
Tfl\x, y, fJ. R /Q 2 , Mf/Q 2 ) from (||) becomes nf (i.e. (3q) independent. 

A glance at Eq. (0), where the coupling constant os(/£r) appears under the 
integral sign, reveals that any of the choices of fi R given by (Qnl) leads immediately 
to the problem if the usual one-loop formula for the effective QCD running coupling 
constant is employed. To circumvent this, one can introduce a cutoff in one-loop 
formula with the aim of preventing the effective coupling from becoming infinite for 
small gluon momenta. There are number of proposals for the form of the coupling 
constant as(n R ) for small fi R [[i~7|, Q, but its implementation in this calculation 
demands the more refined treatment (see the discussion in j^]). Alternatively, one 
can choose (j. R to be an effective constant by taking /j, R = \H R /- Hence, in this 
work we have replaced the expressions (0j^) by their respective averages 



Owing to the fact that 4> a s{x^ 2 F ) (as well as any other pion DA) is centered 
around the value x = 0.5, we take the average value of the momentum fraction to 
be (x) = 0.5. 

We have shown in || that the results depend very weakly on the choice of 
the factorization scale [If- Actually, taking /i F to be an effective constant, i.e., 
H F = (/J- 2 ?), the only [i F dependence of the results obtained using 4> as (x, /i F ) 
distribution comes from the NLO evolution of the DA which is negligible. In the 
following we take n 2 F = Q 2 . 

4- Numerical results 

We take that a perturbative prediction for F 7r (Q 2 ) can be considered reliable pro- 
vided the corrections to the LO prediction are reasonably small (< 30%) and 
the expansion parameter (effective QCD coupling constant) is acceptably small 
( a s(n R ) < 0.3 or 0.5). The consistency with the experimental data is not of much 
use here since reliable experimental data for the pion form factor exist for Q 2 < 4 
GeV 2 |]l9) i.e., outside the region in which the perturbative treatment based on 
Eq. (Q) is justified. It should also be mentioned that one can find controversial ar- 
guments in the literature [^0] regarding the reliability of the analysis of Ref . jl9| . 




(x) 2 Q 2 , 

<-> 3/2 Q 2 , 

e- 5 / 3 (x) 2 Q 2 . 



(10) 

(11) 
(12) 
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Figure 1: (a) Leading- twist NLO QCD results for Q 2 F 7T (Q 2 ) obtained using the 

(f'asix, n F ) distribution amplitude and the choices of n R given byEqs. 

fl2| ), while /ip — Q 2 and (x ) = 0.5. The experimental data are taken from jnj. (b) 

The ratio fP (Q 2 ) / Fi 0) (Q 2 ) obtained using the same DA and the same choices 

for and y? F scales as in Fig. la. The shaded area denotes the region of 

predictions which corresponds to \F^ (Q 2 ) j F^ {Q 2 )\ < 30%, and as(n R ) < 0.5 

(a s (/4)<0.3). 



The new data in this energy region are expected from the CEBAF experiment 
E-93-021. 

Numerical results of our complete NLO QCD calculation for the pion form 
factor, Fk(Q 2 ), obtained using the <j) a8 {x,n 2 F ) distribution, with (i F — Q_ and 
different choices for the renormalization scale /i R given by (||) and (|Io|-|l^), are 
displayed in Fig. |l|a (in our calculation we take ^-jjg = 0.2). The ratio of the NLO 

to the LO contribution to F V (Q 2 ), i.e., F^ (Q 2 )/ F^ ] (Q 2 ), as a useful measure of 
the importance of the NLO corrections, is plotted as a function of Q 2 in Fig. [l]b. 

The solid curve in Figs. |l|a and b corresponds to the often encountered choice 
^ 2 r = Q 2 - The total perturbative prediction F 7T (Q 2 ) is somewhat below the trend 
indicated by the presently available experimental data. To make a meaningful 
comparison between theory and experiment, reliable experimental data are needed, 
as well as a reliable estimation and inclusion of the soft contributions in a moderate 
energy region. What alarms us is the fact that, the ratio (Q 2 ) I Fi \Q 2 ) is 
rather high and F^^Q 2 ) / F^ 0) (Q 2 ) < 30% is not reached until Q 2 » 500 GeV 2 . 
This result seems to question the applicability of the pQCD to exclusive processes. 
The answer to this problem lies in the previously stated inappropriateness of the 
choice = Q 2 . Namely, owing to the partitioning of the overall momentum 
transfer Q 2 among the particles in the parton subprocess, the essential virtualities 
of the particles are smaller than Q 2 , so that the "physical" renormalization scale, 
better suited for the process of interest, is inevitably lower than Q 2 . 
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The results displayed in Fig. |lja show that the total prediction for the pion form 
factor Q 2 F n (Q 2 , fj, 2 ^, fi F ) depends weakly on the choice of n 2 R . This is a reflection of 
the stabilizing effect that the inclusion of the NLO corrections has on the LO pre- 
dictions. By contrast, the results for the ratio F^\Q 2 , y? F ) j F^\Q 2 , n 2 R , ijl 2 f ) 
are sensitive to the choice of From the results displayed in Fig. |l|b we find 
that by choosing the renormalization scale related to the average virtuality of 
the particles in the parton subprocess or given by the BLM scale, the size of the 
NLO corrections is significantly reduced and reliable predictions are obtained at 
considerably lower values of Q 2 , Q 2 < 100 GeV 2 . 



5. Conclusions 

We conclude by stating that, regarding the size of the radiative corrections, the 
sHSA can be consistently applied to the calculation of the pion electromagnetic 



form factor. Further investigation of the scale fixing problem 21 as well as the 
calculation of the NLO prediction in the mHS A (the difference should be important 
only in the region of a few GeV) remain challenges to future work. 
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